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Abstract. Neutrinoless Double Beta Decay (0^/3/3) is presently the only known experiment to 
distinguisch between Dirac neutrinos, different from their antiparticles, and Majorana neutrinos, 
identical with their antiparticles. In addition Ou/3/3 allows to determine the absolute scale of the 
neutrino masses. This is not possible with neutrino oscillations. To determine the neutrino masses 
one must assume, that the light Majorana neutrino exchange is the leading mechanism for Qvj3j5 
and that the matrix element of this transition can ba calculated reliably. The experimental 0vj3j3 
transition amplitude in this mechanism is a product of the light left handed effective Majorana 
neutrino mass and of this transition matrix element. The different methods, Quasi-particle Random 
Phase Approximation (QRPA), Shell Model (SM), Projected Hartree-Fock-Bogoliubov (PHFB) and 
Interacting Boson Model (IBM2) used in the literature and the reliability of the matrix elements in 
these approaches are reviewed. In the second part it is investigated how one can determine the leading 
mechanism or mechanisms from the data of the decay in different nuclei. Explicite expressions 

are given for the transition matrix elements, is shown, that possible interference terms allow to test 
CP (Charge and Parity conjugation) violation. 



PACS numbers: 21.60.-n, 21.60.Jz, 23.40.-s, 23.40.Hc, 
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1. Introduction 

The neutrinoless double beta decay (Ov/3/3) allows as only known experiment to distinguish between 
Dirac (neutrinos different from the antiparticles) and Majorana neutrinos (neutrinos identical with 
its antiparticles apart of a possible phase). It also provides a method to determine the absolute scale 
of of all three neutrino masses in connection with neutrino oscillation data. Neutrino oscillations 
give only the differences of the masses squared. The determination of the masses is possible, if one 
assumes, that the light left handed Majorana neutrino exchange is the leading mechanism for the 
neutrinoless double beta decay and one is also able to calculate reliably the transition matrix element. 
The experimental amplitude called normally T 0u is the product of the effective Majorana neutrino 
mass and a transition matrix element Ai® u . To determine the absolute masses the matrix element 
is as important as the data for the Ov(3f3 transition. The different methods used to calculate 
these matrix elements are presented and compared with their advantages and their drawbacks [lj. 
The methods are the Quasi-particle Random Phase Approximation (QRPA) for spherical (chapter 
2) and deformed (chapter 3) nuclei [H EJ El EJ EJ E] , the She11 Model ( SM ) [3 El E], the Projected 
Hartree-Fock-Bogoliubov (PHFB) approach [TUJ [TTJ [TJl [131 E] an d the Interacting Boson Model 
(IBM2) [15] . (In ref. [13] Pradfully Rath et al. corrected a missing factor 2 and thus in all older 
publications the matrix element have to be multiplied in their work [HI [12] by a factor 2.) 
In the second part the assumption, that the light left handed Majorana neutrino exchange is the 
leading mechanism is not assumed and possibilities to determine the leading mechanism are given 
[Toj [T71 [18] even in cases, where two equally strong mechanisms interfere. There the relative phase 
can test the CP (combined Charge conjugation and Parity) conservation or violation due to relative 
Majorana phases. For CP conservation the strength coefficients r\ must be real and thus the relative 
phase angle can only be zero or 180 degrees. 

Figure [H shows the diagram for the neutrinoless double beta decay of \%Ge^ to USe^ through the 
intermediate nucleus I^As^- 

In addition to the light left handed Majorana neutrino exchange, one has other possible 
mechanisms as cause for the neutrinoless double beta decay: Grand Unification (GUT), 
Supersymmerty (SUSY) and extensions to extra dimensions. We shall discuss here extension to 
GUT's and SUSY. 

We consider first a left-right symmetric model of GUT. 

Wi = cos -&GUT + sin -&GUT W 2 = - sin $ GUT + cos $gut (1) 

Here W\ is the usual vector bosons of Rubbia and coworkers of 80.4 GeV mainly responsible 
for the left handed weak interaction. This allows at each vertex in figure [1] left and right handed 
interactions and one handedness at the leptonic, electron-neutrino side and an other handedness at 
the hadronic, neutron-proton side inside each four-vertex. In addition one has SUSY contributions, 
where mainly the threelinear terms are responsible for the lepton number violation [161 El EEE1 [19] . 
This will be discussed in chapter four of this contribution. 

Fermis Golden Rule of second order time dependent perturbation theory yields: 

r ou = f dF (jA ST < f\H w \k >< k\H w \i > 

J k[ 1 4- £ +( 76 Ge) - [Efc(er) + E k {y) + E k (™As)] 1 1 
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Figure 1. The neutrinoless double beta decay of T^Ge^ to 34^642 through the intermediate nucleus 
33^4s43. The neutrino is emitted at the left vertex as an antineutrino v c with positive helicity due 
to the left handed interaction and must be absorbed at the right vertex as a neutrino with negative 
helicity. This is only possible for a massive Majorana neutrino, which violates helicity conservation 
and is a 50% mixture of neutrino and antineutrino. The two emitted electrons violate lepton number 
conservation. This is possible due to the lepton and antilepton mixture of the Majorana neutrino. 
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is often called the gold plated term. If one wants to determine the effective light left handed 
Majorana neutrino mass from the neutrinoless double beta decay, one assumes that this is the leading 
term and one can neglect the rest. But one needs apart of the data a reliable value for the matrix 
element. In the second chapter we will compare the different methods to calculate these matrix 
elements. 



2. The different Many Body Approaches for the Ou/3/3 Matrix Elements. 



The groups in Tubingen, Bratislava and Jyvaskyla [TJ [2j [3j HJ El |6] are using mainly the Quasiparticle 
Random Phase Approximation (QRPA), while the Strasbourg-Madrid group [TJIH1E] uses the Shell 
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Model (SM), Tomoda, Faessler, Schmid and Gruemmer [TD] and Rath and coworkers [TTJ[T21[T3] use 
the angular momentum projected Hartree-Fock-Bogoliubov method (HFB) (See the correction for 
the missing factor 2 in the work of Pradfulla Rath et al. in [13].), and Barea and Iachello [15] use 
the Interacting Boson Model (IBM2, which distinguishes between protons and neutrons). 

The QRPA [H El El H] has the advantage to allow to use a large single-particle basis. Thus, one 
is able to include to each single nucleon state in the QRPA model space also the spin-orbit partner, 
which guarantees that the Ikeda sum rule [20] is fulfilled. This is essential to describe correctly the 
Gamow- Teller strength. The SM [7] is presently still restricted to a nuclear basis of four to five single- 
particle levels for the description of the neutrinoless double beta decay. Therefore, not all spin-orbit 
partners can be included and, as a result, the Ikeda sum rule is violated by 34 to 50% depending on 
the single particle basis used. On the other side, the shell model can in principle take into account 
all many-body configurations in a given single-particle basis. The excited states in the QRPA for 
spherical even- mass nuclei include 'seniority' (the number of broken quasiparticles) states two, six, 
ten, . . . and for the ground state correlations 'seniority' zero, four, eight, . . . configurations. The SM 
takes for the ground state seniority zero, four, six, eight, . . . , and for the excited states seniority 
two, four, six, ...into account. But the numerical results of the shell model show [HE], that in 
agreement with the philosophy of RPA the contributions of seniority 6 configurations are small and 
can be neglected pQ. In QRPA one starts from the transformation to Bogoliubov quasiparticles : 

a\ = Uic\ - ViCj. (6) 

The creation c\ and annihilation operators of time reversed single-particle states cj are usually 
defined with respect to oscillator wave functions [2]. The single-particle energies are calculated with 
a Woods Saxon potential [2]. The single-particle basis can include up to 23 nucleon levels (all single- 
particle states without a core up to the ii 3 / 2 level and even much larger, if deformation is allowed 
[3,[6]) for the protons and also for the neutrons. But the QRPA results for the Ovfifi matrix elements 
turn out to be stable as a function of the basis size already for smaller basis sets (from 6 or 7 levels 
and larger, respectively) in lighter systems. 

The excited states |m) with angular momentum J in the intermediate odd-odd mass nucleus are 
created from the correlated initial and final + ground states by a proton-neutron phonon creation 
operator: 

\m)=Ql\0 + ); HQl\0+) = E m Ql\0+). (7) 

Ql = J2i x « A «- Y ™ A ^ ( 8 ) 

a 

which is defined as a linear superposition of creation operators of proton-neutron quasiparticle pairs: 

4 = [aWkhM, (9) 

For the present presentation the complication of angular momentum coupling, which must and 
is included in the quantitative calculations, is not shown. 

The inverse 0^/3/3 lifetime for the light Majorana neutrino exchange mechanism is given as the 
product of three factors, 

(T^)' 1 = G 0v \M° U U \ 2 -<m v > 2 (10) 
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where G Qu is a calculable phase space factor, M.^ is the Ou/3/3 nuclear matrix element, and < m v > 
is the (nucleus- independent) "effective Majorana neutrino mass" (j3J). 

The expressions for the matrix elements M^f and the corresponding Ou/3/3 transition operators 
are given, e.g., in Ref. [2]: 

M^ v) = M°Sr ~ (—) 2 M° F U - M° T U (11) 
9a 

The SM approach has been applied by the Strasbourg-Madrid group [3] to neutrinoless double 
beta decay [7] using the closure relation with an averaged energy denominator. In this way one 
does not need to calculate the states in the odd-odd intermediate nuclei. The quality of the results 
depends then on the description of the + ground states in the initial and final nuclei of the double 
beta decay system, e.g. 76 Ge — ?- 76 Se, on the nucleon-nucleon interaction matrix elements fitted by the 
Oslo group in neighbouring muclei and on the average energy denominator chosen (fitted) for closure. 
The 0z//3/3 transition matrix element (fTTT) simplifies as shown in equations (5) to (11) of Ref. [8]. Since 
the number of many body configurations is increasing drastically with the single-particle basis, one 
is forced to restrict for mass numbers A = 76 and A = 82 in the SM the single-particle basis to 
1^3/2,0/5/2, lpi/2 and 0g9/2. In ref. [8] the 82 Se nucleus is calculated in the SM for five basis single- 
particle levels including also 0/7/2- For the mass region around A = 130 the SM basis is restricted 
to 0(77/2, 1(^3/2, 1^5/2, 2si/ 2 and O/tu/2 levels. The problem with these small basis sets is that the 
spin-orbit partners O/7/2 and Ogj/2 have to be omitted [5|. The SM results then automatically violate 
the Ikeda Sum Rule (ISR) [20J, while the QRPA satisfies it exactly. The Ikeda sum rule is: 



S_- S + = 3(N - Z), (12) 

A A 

s - = E<°* + iE(-)" o "-/»(*)*+(*)]E ff M(o*-(0]|of>, (13) 

/1 k I 



For S + the subscripts at the isospin rising and lowering operators are exchanged. 

Figure [2] shows the QRPA contributions of different angular momenta of the neutron pairs, which 
are changed in proton pairs with the same angular momenta. In figure [2] the left bar is the result 
for 82 Se obtained with the single-particle basis 1^3/2,0/5/2, lpi/2 and 0(79/2 used in the SM. The ISR 
is exhausted by 50% only. The second bar from the left represents the result with addition of the 
I/7/2 level. The ISR is exhausted by 66%. The third bar from the left shows the result obtained by 
inclusion of both spin-orbit partners O/7/2 and 0(7 7 / 2 missing in the four level basis of the SM. The 
ISR is 100% fulfilled. For the right bar the basis is increased to 9 single-particle levels for neutrons 
and protons (0/ 7/2 , lp 3 /2, O/5/2, lpi/2, 0# 9 /2, W5/2, 2si/ 2 , l(i 3 /2, 0^7/2)- 

In the last ten years P. K. Rath and coworkers [TH [121 EE] have published a whole series of 
papers (see references in Ref. [12]) on 2v(3(3 decay and, since 2008, also on Qv/3/3 decay, in which they 
used a simple pairing plus quadrupole many body Hamiltonian of the Kumar and Baranger type [21] 
to calculate the neutrinoless double beta decay transition matrix elements with angular momentum 
projection from a Hartree-Fock-Bogoliubov (HFB) wave function after variation. Schmid [22] did 
show, that with the assumption of a real Bogoliubov transformation (real coefficient A and B), axial 
symmetry 

M 
i=i 
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Figure 2. (Color online) Contributions of the transforming neutron pairs with different angular 
momenta J T to the total Mj" calculated within the QRPA and different basis sizes for the 0^/3/3 
decay 82 Se— >- 82 Kr. The left bar is calculated with the same basis of four levels, lps/2, 0/5/2, lpi/2 
and Ogg/2, used in the shell model calculations [7_l[8j[9]. The Ikeda Sum Rule (ISR) [20] is exhausted 
by 50%. The second bar from the left includes in addition the I/7/2 level, one of the two missing 
spin-orbit partners given for the 82 Se nucleus in ref. [8] for the shell model. The ISR is exhausted by 
66%. The third bar from the left includes both missing spin-orbit partners O/7/2 and Ogi/2 amounting 
in total to 6 single-particle levels. The ISR is fulfilled by 100%. This leads to the increase in the 
neutrinoless matrix element from 1.12 to 4.07. The right bar represents the QRPA result with 9 
single-particle levels (1/7/2, 2p 3/2 , 1/5/2, %Pi/2, lffg/2, 2d 5/2 , 3s 1/2 , 2d 3/2 , lflY/2-)- The matrix element 
gets only slightly increased to 4.27. The spin-orbit partners are essential to fulfill the Ikeda Sum Rule 
(ISR). In all four QRPA calculations the QRPA "renormalization" factor g pp (given in the figure) of 
the particle-particle strength of the Bonn CD nucleon-nucleon interaction is adjusted to reproduce 
the experimental 2^/3/3 decay rates. 



and no parity mixing, only + , 2 + , 4 + , . . . . nucleon pairs and excited states are allowed (See eqn. 
(4.2.3) on page 603 of ref. [22]). Rodriguez and Martinez-Pinedo start with the projected HFB 
approach but allow admixtures of different deformations using the Generator Coordinate Method 
(GCM) and the Gogny force [23]. The QRPA and the SM do not have this restriction like PHFB 
and also in the PHFB with the deformation GCM extension (GCM-PNAMP) [TJ]. 
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Figure 3. Contributions of neutron pairs with different angular momenta to the neutrinoless double 
beta decay transition matrix elements for 76 Ge — > 76 Se calculated from a HFB wave function with 
angular momentum and proton and neutron particle number projection before variation. The Fermi, 
the Gamow- Teller and the total contribution including the tensor part as defined in eq. (|lf \ arc 
separately given. The nucleon-nucleon interaction is an improved Gogny type force |23j . The 
deformations ftce = —0.08 and f3s e = 0.11 correspond to the minima of the projected HFB total 
energy. The results for the transition matrix elements are qualitatively and almost quantitatively the 
same for the experimental deformation from the Coulomb reorientation effect: j3o e = 0.16, Pse = 0.10 
and also for different forces. The angular momenta of the neutron pairs are in the PHFB approach 
with axial symmetry, real coefficients and no parity mixing restricted to + , 2 + , 4 + , . . .. In addition the 
contributions of higher angular momentum neutron pairs 2 + , 4 + , . . . are drastically reduced compared 
to the QRPA and the SM. 



Figure |3] shows on the other side, that the projected HFB approach is restricted to contributions 
of neutron pairs with angular momenta + ,2 + ,4 + , .... In addition, one sees that the contributions 
of transition of higher angular momentum neutron to proton pairs 2 + , 4 + , . . . are drastically reduced 
compared to the QRPA and the SM see fig. [3j The reason for this is obvious: in a spherical nucleus 
the HFB solution contains only seniority zero and no stronger higher angular momentum pairs. The 
double beta decay system ^Ge^ — >H Se^ has only small deformations and thus a projected HFB 
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state is not able to describe an appreciable admixture of higher angular momentum pairs for + — > + 
transitions as can be seen in ref. [22]. The higher angular momentum contributions increase drasti- 
cally with increasing intrinsic quadrupole and hexadecapole deformations of the HFB solution. 

The results in figure [3] are calculated in.Tuebingen by K.W. Schmid [22] within the HFB 
with angular momentum and particle number projection before variation with an improved Gogny 
force [23] adjusted in a global fit to properties of many nuclei. 

To have also 1 + , 3 + , 5 + , . . . neutron pairs contributing one has to use a Bogoliubov transformation 
with complex coefficients A and B (fl4|) . To have also 0~, 1~, 2~, 3~, 4~, 5~, . . . one has to allow 
parity mixing in the Bogoliubov transformation ( 1T4"1) . But even allowing all different types of angular 
momentum and parity pairs one would still have an unnatural suppression of the higher angular 
momenta especially for smaller deformations. This handicap could probably be overcome by a multi- 
configuration HFB wave function [22] with complex coefficients and parity mixing in the Bogoliubov 
( fH|) transformation. 

The IBM (Interacting Boson Model) [T5j can only change + (S) and 2 + (D) fermionic pairs 
from two neutrons into two protons. In the bosonization to higher orders this leads to the creation 
and annihilation of up to three "s" and "d" boson annihilation and creation operators in Ref. [13] . 
But all these terms of equation (18) of reference [13] originate from the annihilation of a + (S) or a 
2 + (D) neutron pair into a corresponding proton pair with the same angular momentum. The higher 
boson terms try only to fulfill the Fermi commutation relations of the original nucleon pairs up to 
third order. The IBM can therefore change only a + or a 2 + neutron pair into a corresponding 
proton pair. 

3. Including the Nuclear Deformation in QRPA. 

Table 1. The values of the deformation parameter of the Woods-Saxon mean field /?2 for initial (final) 
nuclei fitted to reproduce the experimental quadrupole moment (labeled as "1"). The spherical limit is 
labeled as "0" . The particle-particle strength parameters g pp are listed. They multiply the Briickcner 
particle-particle nucleon-nucleon G matrix elements of the BONN CD force. They are fitted to the 
2i//3/3 decay half lives. The axial charge is assumed to be the vacuum value qa — 1.25. The particle- 
hole strength parameter, with which the nucleon-nucleon BONN CD Bruckner particle-hole matrix 
elements are multiplied, is adjusted to the excitation energy of the Gamow- Teller resonance in the 
intermediate nucleus g p h = 0.90. The BCS overlap factor (BCSf\BCSi) between the initial and final 
BCS vacua is given in the last column. 



initial (final) 






Qpp 


(BCS t \BSC f ) 


nucleus 










76 Ge ( 76 Se) 


0.10 (0.16) 




0.71 


0.74 




0.0 (0.0) 


"0" 


0.68 


0.81 


150 Nd (150g m ) 


0.240 (0.153) 




1.05 


0.52 




0.0 (0.0) 


"0" 


1.01 


0.85 


160 Gd (160 Dy ) 


0.303 (0.292) 




1.00 


0.74 
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Table 2. The total calculated nuclear matrix elements (NME) Ml" for Ov/3/3 decays 7b Ge-> 76 Se 
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,150c 
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Dy. The BCS overlaps from table 1 are taken into account. In the last 



two columns the 0vf3f3 matrix element M^ 1 and the half-lives for assumed < m v >= 50 meV are 
shown. 
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T°; 2 ■ 

< m v >=50 meV 


76 


"0" 


1.25 
1.25 


4.69 
5.30 


7.15 
5.60 


150 


"0" 


1.25 
1.25 


3.34 
6.12 


0.41 
0.12 


160 




1.25 


3.76 


2.26 



We have also calculated [5J E] the transition matrix elements of the light left handed Majorana 
neutrino exchange Ai® u . Different deformations for the inital and the final nuclei are allowed. The 
BCS overlaps are listed in table 1. The quadrupole deformations are taken from the reorientation 
Coulomb excitation of the 2 + states. 

Q reorientation 



02 



(15) 



5 Z <C T ^charge 

The deformation reduces the matrix elements [5], [6] in 76 Ge slightly from 5.30 to 4.69 by 10% 
only. This is within the error of the matrix elements (see figure E]). But the reduction of the matrix 
elements is severe in strongly deformed systems with different deformations for the initial and the 
final nuclei. In the system 150 Nd — > 150 Sm the matrix element is reduced from 6.12 to 3.34 (see table 
2) and in the strongly deformed system im Gd — > mo Dy (see table 1) one obtains a matrix element 
of 3.76 (see table 2). The single nucleon basis in these deformed calculations are determined in a 
deformed Woods-Saxon potential. The results are then expanded into a deformed oscillator basis 
with the same deformation parameter and the appropriate oscillator length in seven oscillator shells 
[5], [6]. The deformed result for 150 Nd — > 150 Sm is included in figure [51 



4. How to find the Leading Mechanisms for the Neutrinoless Double Beta Decay? 

Normally one assumes, that the first term of eq. ([3]) is the leading one and with the experimental data 
and the matrix element for the light left handed Majorana neutrino exchange M°J one can determine 
the effective Majorana neutrino mass (jlj). But in Grand Unification (GUT) and Supersymmetry 
(SUSY) additional mechanisms for the neutrinoless Double Beta Decay (Oz^/3/3) are possible. 
The matrix element of this gold plated term is proportional to: 

{light PL) OC U ek ■ ^J^^ ■ U ek = J 2 £ " 1^1 W (16) 

fc=i,2,3 n kv y fc=12)3 

The exchange of a heavy left handed Majorana neutrino: 

N e = £ U»N k * e^WR-Nk (17) 

fe=l,...6 fc=4,5,6 
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Figure 4. Diagram for the light Majorana neutrino exchange. The two neutrons in the initial 
nucleus in the ground state + , which change into two protons are are coming from the left and are 
characterized by the two down quarks d in the two neutrons, which change into two up quarks u in 
two protons in the final nucleus. The vector bosons Wl mediating the left handed weak interaction 
are coupled by the expansion coefficients U e k ( [5J to the neutrino mass eigenstates rrikv ■ 



with a k the Majorana phases for these heavy left handed Majorana neutrinos. 



{heavy N L ) cc £ U% ■ P L 1 P L ■ 

i a? a a M kN 



u. 



N 
ek 



fc=4,5,6 



E ^ ■ \ U e k \ 2 /M kN (18) 

fc=4,5,6 



The lepton number and R-parity violating contributions in SUSY are the trilinear terms and 
the coupling of the lepton superfields to the Higgs particle. 

W$ = X ijk ■ Li ■ Lj ■ E c k + \' ijk ■ U ■ Qj -Dl + fMi-Li-Hi (19) 

The lepton L, E and quark Q, D left (L) and right (R) handed superfields are defined as: 



e 
v 



kL 



Eh 



kR 



(U\ 

d 

u 



D k 



(20) 



kR 



kL 



The indices i,j,k run over the three families for leptons: e, li, t and for quarks: d, s, b. The subscripts 
L and R characterize left and right handed superfields. The tilde indicates SUSY particles like 
selectrons , sneutrinos and squarks. 
The inverse half life is given by: 



In 2 



1( ,0i/ 



VnlM% l + Vx'M^Y + \vnr\ 2 \M% r \ 2 \ (21) 
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Figure 5. (Color online) Neutrinoless double beta decay transition matrix elements for the different 
approaches: QRPA [3] , the SM El [8] , the projected HFB method [13], the projected HFB with 
the Generator Coordinate Method (GCM) with deformations [14] (GCM-PNAMP) and the IBM2 [15]. 
The error bars of the filled circles for the QRPA are calculated as the highest and the lowest values 
for three different single-particle basis sets, two forces (Bonn CD and Argonne V18) two different 
axial charges gA — 1-25 and the quenched value gA — 1-00 and two different treatments of short 
range correlations (Jastrow-like [25] and the Unitary Correlator Operator Method (UCOM) 26 ). 
The radius parameter is as in this whole work ro = 1.2 fm. The triangle with the tip up are the 
SM results [7, 9, 8 . The triangle with the tip down represent the transition matrix element of the 
Interacting Boson Model 2 (IBM2) [15]. The squares have been calculated by Pradfulla Rath and 
coworkers [13] with the correction of the factor 2 from December 2010 included, with which all previous 
results of Rath et al. have to be multiplied [UJ [12l [J3] . The star (GCM-PNAMP) is a projected HFB 
calculation with the Gogny force [53] by Rodriguez and Martinez-Pinedo [T3] allowing for different 
deformations with the Generator Coordinate Method (GCM). 
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(23) 
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otik = oti - a k 



(24) 



We restrict here for SUSY to the trilinear terms ( fT9l) .They are lepton number and R parity 
violating. They can contribute by gluino or by neutralino exchange. The strength parameter n\> for 
gluino exchange contains the coupling constant X' 211 the gluino mass m g and the SUSY left and right 
handed up and down squark masses. 



For neutralino exchange one obtains a corresponding expression [161 13 US]. Both exchanges 
can be summarized under the parameter ny for the phenomenological analysis. 

To test, if the light left handed Majorana neutrino exchange is the leading mechnism, one 
meeds at least experimental data of the neutrinoless double beta decay and reliable transition matrix 
elements in two systems. The light Majorana neutrino exchange is represented by the gold plated 
term, which is the first on the right hand side of eq. ([3]) . If the measurements and the matrix element 
are reliable enough [18] and the light Majorana neutrino exchange is indeed the leading mechanism, 
one should in both and also in all other systems obtain the same effective Majorana neutrino mass. 
If two mechanisms are at the same time contributing, one must distinguish between non-interfering 
and between interfering mechanisms. The light left handed Majorana neutrino and the heavy right 
handed neutrino exchange have negligible interference [T7] . 



To determine the absolute values of the two strength parameters r\ v and n^R one needs at least 
two decay systems i. To verify, that these are indeed the leading mechanisms one needs at least 
a measurement in one additional decay system i. But if one forms ratios of half lives using the 
Tuebingen matrix elements for the light Majorana neutrino exchange and the heavy right handed 
neutrino exchange one otaines a very restricted allowed interval for these ratios [T7] . 



The dependence of this ratios on the different parameters of the nuclear structure calculation is 
very minor. Due to the ratios the dependence on most changes drop approximately out. The ratios 
( |27l) are calculated for the axial charge g& = 1.25. But the quenching of this value to qa = 1-00 has 
only a minor effect |17j . 

If the two leading mechanisms like the light Majorana neutrino exchange and the SUSY mechanism 




(25) 




(26) 




(27) 
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with gluino or neutralino exchange can interfer, the situation is a bit more complicated: Let us 
assume the relative phase angle of the complex strength parameters r\ v and rjy is $ Ut \>. The inverse 
of the half life time the phase space factor G® u , 2 (Eo, Z) is then: 

T o» /2 . G l {EojZ) = M 2 (MtT + M 2 (m^) 2 + 

cos0 v ,x>'\r }v \-\r ]x ,\-M%-M%' x r, (28) 

One needs three decay systems to determine the absolute values of the parameters r] v , r/y and 
the relative phase angle fi^y. At least one additional system is needed to verify, that indeed these 
two mechanisms are the leading ones. Again the ratios of the half lives are allowed to lie only in 
narrow regions [17]. If this is not the case, the chosen mechanisms are not the leading ones [T7] . 
With CP conservation the strength parameters r] must be real and thus the relative phase angle is 
zero or 180 degrees. So the determination of "& Ut \i allows to test CP conservation or violation. 

5. The effective Majorana Neutrino Mass. 

Before we summarize the results let us assume Klapdor-Kleingrothaus et al. [21] have indeed 
measured the neutrinoless double beta decay in 76 Ge, although the general belief is, that this still 
needs confirmation. From the half life given by Klapdor et al. [21] one can derive with our matrix 
elements the effective Majorana neutrino mass (jl]). 

Ty 2 ( 76 Ge, exp Klapdor) = (2.23 + 0.44 - 0.31) ■ 10 25 [years]; (29) 

With our matrix elements one obtains the effective light left handed Majorana neutrino mass 
under the assumption, that the light Majorana exchange is the leading mechanism. 



< m u >= 0.24[eV] (exp ± 0.02; theor. ±0.01) [eV] (30) 

The uncertainty (error) from experiment is 0.02 [eV], while the theoretical error originates from 
the uncertainties of the QRPA matrix elements as indicated in figure El The theoretical error is 0.01 
[eV]. 

6. Conclusions 

Let us now summarize the results of this contribution: 

The Shell Model (SM) [7J El [8] is in principle the best method to calculate the nuclear matrix 
elements for the neutrinoless double beta decay. But due to the restricted single-particle basis it 
has a severe handicap. The matrix elements in the 76 Ge region are by a factor 2 smaller than the 
results of the Quasiparticle Random Phase Approximation (QRPA) [TJ [2j [31 H] , the projected Hartree 
Fock Bogoliubov approach [131 E] an d the Interacting Boson Model (IBM2) [15]. With the same 
restricted basis as used by the SM the QRPA obtains roughly the same results as the SM (figure [2]), 
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but the Ikeda sum rule [20] gets strongly violated due to the missing spin-orbit partners in the SM 
single-particle basis. 

The angular momentum projected Hartee-Fock-Bogoliubov (HFB) method [H] is restricted in its 
scope. With a real Bogoliubov transformation without parity mixing and with axial symmetry (fI4|) 
one can only describe neutron pairs with angular momenta and parity + ,2 + ,4 + ,6 + , . . . changing 
into two protons for ground state-to-ground state transitions. The restriction for the Interacting 
Boson Model (IBM) [15] is even more severe: one is restricted to + and 2 + neutron pairs changing 
into two protons. 

A comparison of the 0u(3(3 transition matrix elements calculated recently in the different many 
body methods: QRPA with realistic forces (CD Bonn, Argonne V18), SM with nucleon-nucleon 
matrix elements fitted in neughbouring nuclei, projected HFB [13J with pairing plus quadrupole 
force [21], projected HFB with the deformation as Generator Coordinate (GCM+PNAMP)[14] and 
with the Gogny force [23] and IBM2 [15] is shown in Fig. |5j 
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